Electronic Transport through Magnetic Molecules with Soft Vibrating Modes 
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The low-temperature transport properties of a molecule are studied in the field-effect transitor geometry. 
The molecule has an internal mechanical mode that modulates its electronic levels and renormalizes both the 
interactions and the coupling to the electrodes. For a soft mechanical mode the spin fluctuations in the molecule 
are dominated by the bare couplings while the valence changes are determined by the dressed energies. In this 
case, the transport properties present an anomalous behavior and the Kondo temperature has a weak gate voltage 
dependence. These observations are in agreement with recent experimental data. 

PACS numbers: 72.15.Qm, 73.22.-f 



The recent development of molecular transistors (MT) has 
created new scenarios for the study of coiTelation effects in 
nanoscopic systems. These devices have attracted a lot of in- 
terest due to their potential application in nanoelectronics and 
their rich variety of behavior Molecular transistors consist 
of a small molecule connecting two electrodes and, in most 
cases, a gate electrode is used to control the molecule's charge 
electrostatically. The transport properties of MT show signa- 
tures of strong electronic correlations as Coulomb blockade 
^ and the Kondo effect 121 Isl Jj, Dl] similar to those observed 
in quantum dot devices (QD) OaLZl]. A remarkable difference 
between QD and MT is the coupling, on the latter, of the elec- 
tronic degrees of freedom with a discrete set of mechanical 
modes fflllS. 

For the simplest case of a linear modulation of the 
molecule's electronic levels by a single vibration mode, sev- 
eral effects are predicted to occur The Franck-Condon renor- 
malization of the molecule-electrodes coupling is expected to 
produce a suppression in the sequential and cotunneling trans- 
port through the molecule lIlU Il2l 11311 . The reduction of the 



effective Coulomb repulsion in the molecule may lead to an 
effective e-e attraction il4ll5lll6ll and a strong sensitivity to 
gate voltage. 

In the repulsive e-e interaction regime the spin-Kondo ef- 
fect dominates the low-temperature physics. The Kondo effect 
generates in these devices an increase of the conductance with 
decreasing temperature and a zero-bias peak in the differential 
conductance. These observations are a direct consequence of 
the formation of the Abrikosov-Suhl or Kondo resonance be- 



low the Kondo temperature Tk 01711 . 

In some MTs based on organometallic molecules, an 
anomalous gate voltage dependence of the transport proper- 
ties has been reported 111 S |5|l. In the transition metal com- 
plexes studied by Yu et al. ||5|], Tk depends weakly on the 
applied gate voltage and shows a rapid increase only close 
to the charge degeneracy points. Moreover, the edges of 
the Coulomb blockade diamonds are not well defined in the 
Kondo charge state. Such a behavior is inconsistent with the 
usual theory based on the Anderson model. 

In this Rapid Communication we present a study of the 
Anderson-Holstein model showing results obtained with the 



numerical renormalization group (NRG) OlSl Il9ll . We find 
that, as the frequency of the vibrating mode decreases, an 
anomalous gate dependence of Tjc and of the transport prop- 
erties emerges. This effect arises because the soft vibrating 
modes in the MT drive the system into a new regime where 
the characteristic energy scales for spin and charge fluctua- 
tions are not related as in the conventional theory of the Kondo 
effect. 

The model Hamiltonian is H = Hm + He + Hme where 
the first two terms describe the isolated molecule and the elec- 
trodes, respectively, and the last term describes their coupling. 
We have 

Hm = ednd + Undi^ndi-\{nd~l) {a+a^)+LOoa^ a (1) 

He ^ ^ £ak clkaCaka^ (2) 



k,(T,a 



Hme = ^ Vak{dl Caka + cl^„d„ 



(3) 



k.(T,a 



Here Ud = Udi + Udi, Uda ~ d^d^r, d\ creates an electron at 
the molecular orbital with energy Ed, U is the intramolecular 
Coulomb repulsion and cj^.^.^ creates an electron in the mode k 
of electrode a = L, R. The operator a^ creates an excitation 
of the vibronic mode with energy ujq. The Fermi energy Ep 
is set to zero, h = 1, and all energies are in units of half the 
electrodes' bandwidth. For the sake of simplicity, from here 
on we consider a symmetric molecule with identical L and R 
electrodes (Sak = £k) and take VLk = Vnk = Vhyb- 

The energy spectrum of the isolated molecule (Vhyb — 0) 
can be readily obtained. The energies for the ti^ = 0, 1 and 
2 charge states, with m vibron excitations, are i?o,m = ^^^/ 
cjQ + mijj{), Ei^ra = £d + mujo, and £^2,™ = -A^/cjq + 2£d + 
U + mujo, respectively. The corresponding states are indi- 
cated as |0, m), \a,m), and \2,m) where a is the spin index 
in the Ud ~ 1 charge sector — note that the vibronic states for 
different charge states correspond to different equilibrium po- 
sitions of the coordinate associated with the vibronic motion. 
It is convenient to define the effective single-electron energy 
and e-e repulsion as Sd — -Ei.o ~ Eq^q ~ Sd + A^/wq and 
U = £^2,0 + £^0,0 - 2Ei,o = U ~ 2X'^/uio, respectively. The 
charge degeneracy points of the isolated molecule are given 
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FIG. 1: (color online) Spectral density Ad{uj) (in units of I/ttF). 
Panel a) to f) ^d(cj) vs tj for [/ = 0.2, Vhyt = 0.1/^2, g = 
6.25, A — 0.03 and different values of Sd- g) Color map of Ad{iLj) 
constructed by linear extrapolation of the NRG results obtained for 
different values of Ed (indicated by arrows). Thick dashed lines are 
a guide to the eye of the form Sd = i^ + a and the horizontal line 
indicates the charge degeneracy point for the isolated molecule (sd ~ 
0). Values of Ad{u>) larger than 0.5 are white in the color scale. 



hy Ed — and Sd + U — 0. When the molecule is coupled 
to the electrodes, the charge fluctuations are controlled by the 
Franck-Condon (FC) factors 7„^„i = | (c, n [dj. I 0, m) I . We 
define 7„i = 7o,to ~ e^^g™ /m\, where g — (A/cJo)^ is 
a dimensionless parameter For g ^ 1, transitions between 
low-lying states are exponent ially suppressed and the system 
is in the FC blockade regime illlll2ll . 

We now present numerical results for the molecular spectral 
density Ad{Lo) which determines the transport properties of 
the MT. Ad{io) is given by — (l/7r)ImGdd('^), where Gdd{^) 
is the electronic Green's function of the molecule in the pres- 
ence of the electrodes. The zero-temperature NRG results for 
Ad{iLj) are shown in Fig. [T]for different values of the pa- 
rameters in a regime of strong FC effect {g — 6.25). In the 
electron-hole symmetric case [Sd — —U/2, Fig. [T^)], Adiio) 
shows broad structures at the bare energies Sd and Sd + U and 
a well defined Kondo peak at the Fermi level. Although not 
well resolved in the figure, peaks are also obtained at Id and 
£d + U — these are the first FC sidebands. The evolution of the 
spectral density Ad(aj) as the molecular energy is shifted up- 



wards presents interesting features. First, the broad peaks and 
the first FC peaks-which are better resolved for small | Ed \ [see 
Fig. [I];)]- shift with the molecular energy, while the Kondo 
temperature, as given by the width of the Kondo peak, is not 
very sensitive to it. As Id crosses Ep, the behavior of Ad{uj) 
indicates a rapid change in the occupation of the molecular or- 
bital (ud) ■ There is no crossing of a wide resonance through 
the Fermi level, as it occurs in the g = case. Here instead, 
there is a rapid transfer of spectral weight from w ~ e^ to 
u! ^ Sd + U, in agreement with the atomic limit results Ill9ll . 
The change in the structure of Ad{io) occurs as the Kondo 
peak evolves into a FC peak at small w. This is better seen 
in the color map of the spectral density shown in Fig. [T^). 
This anomalous behavior has important consequences on the 
thermodynamic and transport properties of the system. 

The occupation of the molecular orbital versus the bare en- 
ergy Ed is shown in Fig. |2^) for different values of g. As g in- 
creases, the width of the region where the magnetic configura- 
tion is stable i{nd) — 1) and the crossover widths between the 
different charge states decrease. The former is due to the re- 
duction of the effective e-e repulsion U while the latter, being 
exponential with g, is a manifestation of FC effects |lUll2ll . 
Similar FC effects are obtained in the spinless case 111 ill . How- 
ever, for the Anderson-Holstein model presented here, the in- 
terplay between e-e and e-v interactions leads to a crossover 
between different charge states that is wider than in the spin- 
less case. Also, for the spinfull model with large g the charge 
degeneracy points (with (ud) = 1/2 and 3/2) are not given 
exactly by e^ = and Id + U — as self-energy corrections, 
due to the hybridization, shift the dressed molecular levels. 

The magnetic susceptibility x{T) of the molecule is shown 
in Fig. I2J5) for different values of Cd- In the Kondo regime, 
the low-temperature susceptibility curves collapse into a sin- 
gle curve when the temperature is properly scaled. To cal- 
culate the Kondo temperature Tk and its dependence on 
the molecular energy £d we use Wilson's criterion 12011 : 
Tkx{Tk) / Iji-b'^ — 0.025. The results for Tk are shown in 
Fig. I2};). As Wo decreases (we keep A fixed and therefore g in- 
creases), Tk vs. Id deviates from the usual behavior: it shows 
a weak dependence on Id at the center of the Coulomb block- 
ade region and a fast increase close to the charge degeneracy 
points. 

In terms of the Kondo coupling J the Kondo temperature 
is given by Tk = De~^^f""' where I? is a high energy cutoff 
and po is the bare electronic density of states of the electrodes 
at the Fermi level. While the value of J obtained using sec- 
ond order perturbation theory in Vhyb underestimates the ratio 
Tk/Tk for small values of wq, it is instructive to use this ap- 
proach for a qualitative interpretation of the numerical results. 
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-£d + ujom Ed + U + LOom 



(4) 



For large g, the FC factor 7,, 
with a width of the order of 1 



is peaked at m ^ m* = g 
'm* . As in the summation of 
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FIG. 2: (color online) a) Occupation of the molecular orbital {rid) 
vs Ed for 3 = 0, 2, 4, 6.25 and 9. b) Susceptibility Tx{T)/^b^ vs 
T in a log-log plot for Sd/U increasing from —0.5 to —0.38, and 
g = 6.25. Dashed line indicates Wilson's value. Other parameters 
as in Fig. [T] c) Kondo temperature Tk /T% \s £d for different values 
of g as indicated in the inset. Open squares are the values of Tk 
obtained by fitting the temperature dependence of the conductance 
for p = 6.25 (see text). 



Eq. dU the denominators are slowly varying around m*, the 
Kondo coupling can be approximated as J ~ 2V^^[l/(— ed+ 

m*ujo) + 1/{U + Ed + m*uJo)]. This shows that the rele- 
vant virtual charge fluctuations have a characteristic energy 
given by the bare molecular parameters ej. = £d — oj^m* 
and £d + U = Erf + [/ + ojQTn* corresponding to the broad 
peaks in Fig. [T^). However, as e^ — > (— C/) the m = 
term in Eq. (|4| diverges with an exponentially small prefactor 
e^^, indicating that the perturbation theory breaks down for 
Ed exponentially close to the charge degeneracy points of the 
isolated molecule. In other words, for small cjq and large g, 
the virtual charge fluctuations leading to the Kondo coupling 
are antiadiabatic while the charge instabilities are controlled 
by the dressed energies Ed and Ed + U. 

We stress that, even far from the degeneracy points, the per- 
turbative approach of Eq. (|4|i gives a quantitative estimate of 
Tk/T^ only for large values of cjq ?t U/2. As ujq decreases 
the states \cr,m) with ni > also participate in the ground 
state and should be taken into account ll23n . We are not aware 
of any reliable method to properly handle such a case analyti- 
cally. 

The linear conductance G of the molecular junction at zero 
bias is given by the spectral density Ad{uj) presented above. 
For symmetric electrodes 12\ i25i1 
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FIG. 3: (color online) a) Temperature dependence of the conductance 
for different values of Sd/U, and g = 6.25. The charge degeneracy 
point is at Ed/U ~ —0.38(7. b) Zero-temperature conductance vs 
Ed- Parameters as in Fig. [T]and different values of g as shown in 
the inset, c) Differential conductance G{V) (in units of 2e^//i) as a 
function of the bias and gate voltages (cf. Eq. Q). 



temperature dependence of the conductance for g — 6.25 
and different values of Ed is presented in Fig. [3^). In the 
electron-hole symmetric case (Ed — —U/2) the characteristic 
Kondo behavior is obtained. As Ed increases, the low tempera- 
ture conductance decreases in quantitative agreement with the 
Fermi liquid zero-temperature results (see below). Fore^ > 
the conductance shows a maximum at a temperature T* con- 
sistent with the energy of the first EC side-band. Within the 
Kondo regime our results for the temperature dependent con- 
ductance are well fitted by the phenomenological expression 
G{T) ^ G;^[l + (2i/''-l)(r/TK)2]-''- with s between 0.22 
and 0.25 1^. The value of Tk obtained with this fitting pro- 
cedure is consistent with Wilson's criterion as shown in Fig. 
|2};). It is important to point out that while for g = the 
phenomenological expression correctly describes the univer- 
sal scaling of the conductance even deep into the intermediate 
valence regime ({ud) ~ 0.5), it fails to do so for g ^ 1 al- 
ready for (ud) — 0.75 [Ed — — 0.39C/ for the parameters of 
Eig.EJi). 

The zero-temperature conductance, shown in Fig. ^) as a 
function of Ed for different values of g, was calculated using 
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T^ind) 



(6) 



where f{uj) is the Fermi distribution and F = 2TipoV^ j^. The 



Note that while at high temperatures {T > Tk) the conduc- 
tance is suppressed by the Coulomb interaction and Franck- 
Condon blockade |13], at low temperatures the Kondo effect 
sets in and the conductance is high. 

A salient feature in the differential conductance of systems 
with strong e-v coupling is the anomalous behavior of the 



edges of the Coulomb blockade diamonds. In some MT f2, 
ISfl the diamond edges associated with the Kondo charge state 
((rirf) ~ 1) are much weaker, or even absent, than the ones 
corresponding to the non-Kondo charge state. In the lowest 
order on the bias voltage, the differential conductance is 



G{V) ~ -TTr[AdieV/2) + Adi-eV/2)] . 



(7) 



Although the results obtained with this equation [see Fig. [3};)] 
should only be taken as a rough estimate, they show that the 
qualitative behavior of the differential conductance is very dif- 
ferent from that of the Anderson model. The main reason is 
that, as discussed above, the valence change is not due to a 
broad resonant level that crosses Ep. The spectral weight is 
rather transferred from below to above Ep and only the nar- 
row peak of the first FC band, that has an exponentially re- 
duced weight (oc e^^), crosses Ep- 

In summary, we have studied a model MT where the elec- 
tronic levels of the molecule are modulated by a well de- 
fined internal vibronic mode. We have shown that the cou- 
pling to a soft phonon mode (wq ^ U) in the FC regime 
(g ^ 1) changes qualitatively the behavior of the MT phys- 
ical properties. The electronic interactions in the molecule 
and the coupling to the electrodes are strongly renormalized. 
While spin fluctuations in the molecule are associated to vir- 
tual processes dominated by the bare electronic energies, the 
valence changes occur when the dressed energies cross the 
Fermi level. This leads to an anomalous gate voltage depen- 
dence of the spectral density and Tk- The theory accounts 
for the observed weak dependence of Tk on gate voltage and 
suggests the possible origin of the anomalous behavior of the 
Coulomb blockade diamond edges [5]. We found that the uni- 
versality characteristic of the Kondo phenomena, is lost in this 
case much before the valence change as compared with the 
usual Anderson model. Finally, we showed that the FC sup- 
pression of the low-bias conductance expected for the high 
temperature regime lll2[|l3ll is not present at low temperatures 
due to the Kondo enhancement of the conductance. 
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bers of CONICET (Argentina). 
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